arXiv: 1508.07114v2 [hep-th] 15 Feb 2016 


Oscillation in power spectrum of primordial gravitational wave as 
a signature of higher-order stringy corrections 

Yong Cai^*, Yu-Tong Wang^^, and Yun-Song Piao^’^^ 

^ School of Physics, University of Chinese Academy of Sciences, Beijing 100049, China and 
^ State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, P.O. Box 2735, Beijing 100190, China 

Abstract 

In low-energy effective string theory, a' corrections involve the coupling of the dilaton field to 
Gauss-Bonnet term. We assume that the dilaton potential is fine tuned so that the dilaton field 
may oscillate rapidly for a while around the minimum of its potential but the inflation background 
is not affected. By numerical method, we find that if the dilaton starts to oscillate at the time 
of about ~ 60 e-folds before the end of inflation, a' correction may bring unusual oscillations to 
the inflationary gravitational wave spectrum, which might be measurably imprinted in the CMB 
B-mode polarization. 
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I. INTRODUCTION 


In the recent years, high-precision CMB observations [1] confirmed the predictions of in¬ 
flation. Meanwhile, the searching for the primordial tensor perturbations, i.e., the primordial 
gravitational waves (GWs) [2],[3], which may be imprinted in the CMB B-mode polarization 
spectrum [4], [5], has been still on road, e.g., [6]. The detection of primordial GWs would 
verify general relativity (GR) and strengthen our confidence in inflation. Inflation occurred 
in the early universe with higher scale. Thus, it could be expected that the primordial GWs 
might encode information of the UV-complete quantum gravity (e.g., [7]), for which the 
most promising candidate is string theory. 

The effective action of string theory reduces to GR at low energy. However, it will acquire 
the higher-order curvature corrections with the a' expansion. In heterotic string theory, the 
leading-order correction is [8] 

~ c\Oi e ^Rq^, (1) 

where Rq^ is the Gauss-Bonnet (GB) term. Mg = is the string scale and 0 is the 

dilaton field. In Type II string theory, it is ~ a^R^, e.g., [9], see also Ref. [10] and references 
therein. 

The effects of a' corrections are actually negligible at low energy Mg. Thus the 
common expectation is that these characters of string theory can be hardly directly tested 
in laboratory experiments. However, inflation, which occurred in higher energy scale, is a 
natural laboratory for testing higher energy physics, since the corresponding physics might 
have an influence on the primordial perturbations produced during inflation. The model 
buildings of inflation in string theory has been widely investigated, e.g., [11] and [12],[13] 
for recent reviews. 

The inflationary GWs is a powerful tool to probe the physics beyond GR, since it is 
only affected by the physics relevant to gravity, e.g., the change of the propagating speed 
of GWs [14]. Though the primordial scalar perturbation is affected also by the modified 
gravity, it is mainly affected by the dynamics of infiaton field itself, which will interfere the 
identifying of corresponding signatures. In Refs. [15],[16],[17],[18], the non-Gaussianities of 
inflationary GWs from modified gravity were studied. However, such non-Gaussianities are 
far below the sensitivity of on-road measurements. Gomparatively, the power spectrum of 
primordial GWs is observationlly promising. Thus a significant asking is if a' corrections in 
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string theory may leave measurable imprints in the CMB through their effects on the power 
spectrum of GWs. 

Here, we assume that the dilaton potential is fine tuned so that the dilaton field may 
oscillate rapidly for a while around the minimum of its potential but the inflation background 
is not affected. By numerical method, we find that if the dilaton starts to oscillate at the time 
of about ~ 60 e-folds before the end of inflation, a' correction (1) can lead to the unusual 
oscillations in the power spectrum of primordial GWs. These oscillations may result in 
some obvious wiggles in the GMB B-mode polarization spectrum, which is well within the 
detection of the upcoming experiments. The intensity of the wiggles is determined by the 
string scale Mg and the string coupling g^. 


II. THE SETUP 


The effective action of heterotic string may be [8], 

-S' = ^ y d'^xy/^ + Sinf, (2) 

where = l/M^, and ci = 1/8 is dimensionless. The background is the inflation with 
e = 1, and the field driving inflation is given in Sinf (see, e.g., [13]), which we 

will not involve. The dimensionless scalar 0 is the dilaton field. By variation with respect 
to metric, we obtain the modified Friedmann equations 


a 

a 


6M2 


6^2 ^ ^ 1 ^2 ^ y 2Acia'e-^H^ , 

2 

{Pinf + 3pinf) — + dcicr'e { S ~ ^ ~ 


(3) 


(4) 


6 6 ' \ 2 2H Hcj) 

where pinf and Pinf are the energy density and pressure contributed by infiaton respec¬ 
tively. Here, for simplicity, we set 0 as a ‘non-background’ field, i.e. spectator field, which 
means that it dose not dominate or affect the evolution of the background. This requires 
VMp -C Pinf — and the contribution of cia'e~'^<j)H^ from the GB term to 

must be also negligible, which put the constraint Cia'e~'^<f>H 1. In addition. 


M!gs 


< 1 


(5) 


must be imposed, so that other higher-order a' corrections may be neglected, where gs = e'^ 
is the string coupling. 
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The dilaton 0 is massless in the weak coupling regime, i.e., 0 —1 and gs -C 1. However, 

in the moderately strong coupling regime, 0 ~ — 1 , the potential must exhibit some structure 
to make 0 trapped at the minimum of the potential, so that the equivalence principle is 
preserved at late time [19]. Though the detail of the dilaton potential H(0) is still an issue 
in the study, it is generally thought that V ( 0 ) with a local valley can naturally arise due to 
some nonperturbative effects [ 20 ]. 

Inflation happens when inflaton dominates over other helds. We can see from Eq.(3) that 
it’s unlikely 0 <^ — 1 at the beginning of inflation, since the term with e~'^ will dominate 
over the contribution from inflaton. Thus, it’s natural to assume that 0 locates at the strong 
or moderately strong coupling regime with 0 ~ — 1 at the beginning of inflation, and runs 
towards the regime with 0 < — 1 until trapped by the valley of H(0). In addition, we also 
require that the initial value of H( 0 ) must also be negligible compared with the inflaton 
potential. 

In this setup, the adiabatic scalar perturbation is contributed by inflaton, which is irrele¬ 
vant with the ‘non-background’ field 0. We, following [21], may obtain the quadratic action 
of the GWs mode hij as follows 


^(2) = m2 


dT(fx 


o^Qt 




( 6 ) 


where hij satisfies dihij = 0 and ha = 0 , r is the conformal time, dr = dt/a, a prime denotes 
the derivative with respect to r, ct is the propagating speed of primordial GWs, and QtM^ 
is regarded as effective Planck scale Mpg^j(r), and 


Cp — -|- Ij , Qt — 4:Cia'^H -|- 1, 


(7) 


where ^ = l/e"^, which are consistent with the earlier results in Ref. [22] [23]. Thus, although 
0 dose not dominate the background, due to its coupling to Rq^, h may affect the tensor 
perturbation, see also [17] for similar case but with nonminimal derivative coupling to R. 
Here, the conditions of avoiding the ghost instability are > 0 and Qt > 0. 

Here, if 0 sits in the minimum of its potential or slow rolls, we have Qt/HQt 1 and 
ct/Hct 1, the primordial GWs spectrum will only acquire a small correction [24],[25], 
and its shape will not be altered essentially. However, before resting at the minimum of its 
potential, if 10,0 3> H^, the dilaton held might rapidly roll down towards the minimum of 
its potential and oscillate around it, which will lead to a short-time but drastic variation of 
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the propagating speed of GWs, see (7). We will show that it is this variation, if happened 
during inflation, that will naturally imprint the leading-order string correction (1) to the 
large-scale primordial GWs. 

In the pre-big bang scenario, the a' corrections were applied to solve the singularity 
problem of the big bang model, which will also leave imprints in the primordial GWs [26],[27]. 
However, these imprints appear at small scales, and are not relevant to the CMB polarization, 
unless the pre-big bang phase is followed by an inflationary phase [28],[29], or see [30]. 


III. THE SIGNATURE OF HIGHER-ORDER STRING CORRECTION 


The power spectrum of primordial GWs is calculated as follows. The Fourier series of 
is 

/ ^3 L 

^_^g-»k.x ^ hx{T,k)e\j\k), (8) 

where h\{T,k) = hx{T, k)a\{k) -|- h\(T,—k)a\{—k), the polarization tensors satisfy 

(k) = 0, eSf^(k) = 0, and = 6xx', = eSj^^(-k), the commutation 

relation for the annihilation and creation operators aA(k) and a^(k') is [aA(k),a^, (k')] = 
<^AA'<^*'^Hk — k'). Thus we have the equation of motion for u{t, k) as 

u" + ^ u = 0, (9) 


where 


u{T,k) = hx{T,k) ■ zt, zt 


aMpx/ Qt 
2 


( 10 ) 


Initially, the perturbations are deep inside the horizon, i.e., c^k"^ ^ the initial condition 
is M ~ J The power spectrum of GWs is 


Pt = 


k^ 


A—“h, X 


Ak^ 1 
7r2M2 ■ 


aH/k ^ 1. 


( 11 ) 


Here, since both ct and Qt experience the short-time but drastic variations, analytical 
solution is difficult to obtain. We will first exhibit the numerical results, which help to 
clearly display the nontrivial feature of GWs spectrum induced by a' correction. 
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A. The numerical solution 


To facilitate the numerically solving of the perturbation equation, we introduce the vari¬ 
able a = Ina [31]. Using ^ the evolving equation of 0 and Eq. (9) can be 

rewritten as 


+ 3 + 


H 


+ 24 1 + 


H 


H^cia'e 


+ 14 = 0 

m ’ 


( 12 ) 


^ “^) ““ ^ ~ “]) “ - O' (10) 

where the subscript ‘a’ denotes d/da, and and zt are given by Eqs. (7) and (10), 
respectively. The initial condition of Eq. (13) is written as 


1 1 du 

- 


V 2e"i7 


a<ln(A)’ 


(14) 


where the initial value of ct is set as unity. We assume, without loss of generality, that the 
background of slow-roll inflation is described by constant e, so that the Hubble parameter 
may be written as H{a)/Mp = Ah ■ in which Ah is constant. 

As has been argued, in the moderately strong coupling regime 0 ~ — 1, due to the non- 
perturbative effects, the potential must exhibit some structure to make 0 trapped at the 
minimum of the potential. Initially, 0 might deviate from the minimum of its potential. 
In (2), the dimension of U(0) is the square of mass and 0 is dimensionless. In term of 
Ref. [20], both and are ~ = M/, in which Ve is the volume of extra 

dimensions. Thus with (5) and Qs = <C 1, we have ^ H^. Therefore, 0 will rapidly 
roll down, and oscillate around the minimum of the potential before getting rest. As long 
as ~ ~ M/ are hxed (note that U(0) <C is still reserved), the result is actually 

insensible to other details of U(0). 

Here, for the purpose of numerical simulation, we set the potential locally as U(0) ~ 
, where < —1. To not interrupt inflation, it’s convenient to parameterise the 
potential locally as 


U(0) = Uo • 




(15) 


l + R-M-2Uo(0-</>+)"’ 
where Vq has the dimension of the square of mass, R is a dimensionless constant. The 
purpose of introducing B is to guarantee that V{(/>) doesn’t affect the background evolution 
when 10 — 0*1 is large enough, as we can see from Eq.(3) and (4). We will set B = 10^°, then 
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the value of V{(j)) is smaller than thus is negligible compared with There is 

some hnetunning in the shape of the potential. The shape of dilaton potential is actually 
still an issue in study, and it is generally thought that in heterotic string theory the potential 
with a local valley might arise due to some nonperturbative effects. What we mainly focused 
on is that how the primordial GWs spectrum will be affected, if there is such a potential. 

Because both the contributions from the Gauss-Bonnet and the potential is to move 0 
towards smaller value, 0 won’t be frozen if it located at 10 — 0*1 > 1 initially, and we don’t 
need to create very fine-tuned initial conditions for 0 Below, we will set 0* = —2.6, and 
the initial value of 0 as (j)initiai = —1.5. 

In Fig.l, we plot the evolutions of 0, and Qt- The field 0 rolls down along its potential, 
and gets rest at the minimum of the potential after oscillating around it, since the inflation 
rapidly dilutes its kinetic energy. This evolution induces the oscillations of and Qt- The 
oscillation amplitude of is far larger than that of Qt, since is dominated by 0 while Qt 
is by 0, and 0 can hardly be large enough since -C Pinf — We can see that both 

Ct > 0 and Qt > 0 are always satished in the regime we are interested in, thus there is no 
ghost instability. Here, from the values of relevant parameters, we have H'^/{M^gs) — 0.1, 
which satishes the condition (5). 

It might be concerned that whether the contribution of the spectator scalar held 0 to the 
background evolution can be safely negligible. To verify this issue, we set the contribution 
of 0 in Eq.(3) and (4) as 

= Q 0 ' + G(0) + 24ci«'e-'^i730^ /(16) 


C^2 = 


6 6 




I i 


2H H 4 , 


(17) 


respectively. We plot and with respect to a in Fig. 2, using the same parameters used 
in Fig.l. We can see that the contributions from the terms relevant with 0 are negligible 
compared with that of inhaton in the regime we are interested in. 

In Fig. 3, we plot the ratio of the power spectrum Pt to for different a', Vq and 
e, in which P^^ = ^ 2^2 is that of standard inhation without a' correction. In all cases, 
unusual oscillations in Pt are present, which is obviously a universal character attributed 


^ We plot 0(a) with different initial value (piniuai in Fig. 7, 
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FIG. 1: The evolutions of the dilaton field (j), and Qt for e = 0.003, Ah = 2.72 x 10“®, ci = 1/8, 
a'M2 = 9.6 X 106, Vo/M^ = 6.5 x lO'^, B = lO^o, = -2.6, = -13.4 and cl)initiai = -1-5. 



a 


(b) Cy and Qt 




FIG. 2: The contributions from non-trivial variation of (j) and its derivatives to the background 
evolution, where e = 0.003, Ah = 2.72 x 10“^, ci = 1/8, a'Mp = 9.6 x 10®, Vo/Mp = 6.5 x 10“®, 
B = 10^®, (/* = —2.6, o-initiai = —13.4 and (j)initiai = —1.5. We can see that <C 1 in the regime 
we are interested in. 

to the nontrivial behavior of ct and Qt, or more intrinsically said, the string theory might 
manifest. Though the particle production during inflation may also significantly affect the 
power spectrum of GWs [32],[33], it only leads to a bumplike modification, which is entirely 
different from the behavior of the oscillation showed here. 

In Fig. 4, we plot the CMB lensed BB-mode power spectrum, compared with that of in¬ 
flation without a' correction, in which r = , and both P^^ and P^^ are those of 

standard inflation. Here, we have assumed that the power spectrum of the scalar perturba- 
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FIG. 3: Power spectra of primordial GWs, in which Ah = 2.72 x 10 ci = 1/8, (/>* = —2.6, 
ainitiai = -13.4 and (^initial = -1.5, while e = (0.003, 0.001, 0.003), Vo/M^ = (6.5, 12, 12) x lO'^, 
B = 10^° and a'M^ = (9.6 x 10®, 8.0 x 10®, 9.6 x 10®), for brown solid curve, green dashed curve, 
and magenta dotted curve, respectively. 



Miiltipole I 



Multipole I 


(a) r = 16e = 0.016, Vq/M^ = 6.5 x lO-'^ 


(b) r = 16e = 0.048, Vq/M^ = 1.2 x 10-^ 


FIG. 4: Wiggle features in the CMB lensed B-mode power spectra, in which Ah = 2.72 x 10 ®, 
Cl = 1/8, a'Mp = 9.6 x 10®, B = 10^®, = -2.6, ainuiai = -13.4, (l)iniUai = -1-5, and 

r = . Both P^^ and P^^ are those of slow-roll inflation without a' correction. BKP 

data are from the joint analysis of BICEP2/KeckArray and Planck [6]. 


tion is not affected by the nontrivial behavior of the dilaton 0, and is still ■ In addition, 
we also assnmed that after inflation the dilaton rests at the minimnm of its potential so 
that GR is covered. Thus, the effect of varying ct at late time may be neglected, or see 
[34] [35]. When the dilaton held 0 is hxed, the GB term will not affect the evolution of 
the background, since it is the generalised Euler density of four-dimensional spacetime, and 
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thus is a topological invariant and does not contribute to the equations of motion. 

There are generally a recombination peak at / ~ 80 and a reionization bump at I < 10. 
Both are the imprints of the primordial GWs. The detection of the recombination peak will 
be a confirmation that inflation has ever occurred, while detecting the reionization bump at 
low-/ would help us to understand the physics of pre-inflat ionary universe, e.g.,[36]. When 
the comoving oscillating scale in GWs power spectrum is set at / ~ 80, for r > 0.01 we may 
see obvious wiggles around the recombination peak, which may be tested by the upcoming 
polarization data. 


B. The analytic estimate 


The oscillations in Pt are induced by a' correction, which naturally implies that the 
intensity of the oscillations is closely related to the value of the string scale Mg = Xj^fa'. 
We have observed in Fig. 3 that with larger a', i.e., smaller energy scale of string M^, we 
got larger amplitude of the oscillation in the power spectrum. To acquire more insights on 
relevant physics, we will give an analytic estimate for it. 

To not affect the background evolution, 0^ 12, 1/ <C and AH‘^cia'e~'^(j)a 1 

must be satisfied. Thus we have Qt ~ 1 and ~ 1 — 4:H'^cia'e~‘^(l)aa- In addition, 
we have ~ ^ 1, so when (j) rapidly rolls down and oscillates around the 

minimum of its potential, which corresponds to 0 = 0*, in Eq. (12) dominates (paa, 

i.e., The second term in Eq. (12) is not negligible, which is responsible for 

making the oscillation stop rapidly. Thus our estimate is precise only for the hrst two or 
three oscillations. However, it is enough for us to grasp the physics behind the oscillating 
Pt- When 0 oscillates around 0*, will acquire an oscillation with same frequency. The 
extreme values of the oscillating amplitude of 0% can be estimated as 

Clmaximin) ~ 1 ± dCiCr'e"'^* | |, (18) 


where e is approximately replaced with e Therefore, for Ml 

cia'eM*\Vl > 0.1, which is required for the spectrum showing itself feature. 

We assume that ct experiences an ideally steplike jump from cti to ct 2 , 
Ref. [14] (see, also [37] for the scalar perturbation), we have 

o _ rjinf f{k,ko,x) 


we may have 


and following 


( 19 ) 
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where 


f{k,ko,x) 






( 20 ) 


for k 3> /co, which oscillates between and 1/x'^, and k^ is the wave number of the GWs 
mode crossing horizon when ct jumps, and x = 

Thus we approximately have 


Pt 


n 


inf 


Cti 


-'T min 



min 

max 


( 21 ) 


This reflects the effect of a', which corresponds to Mg, on the amplitude of oscillation in the 
power spectrum. In turn, we have 


m; 


a'Ml 


.. 9s 

'{PrlPT^UaP 

.{PT/PT^)miu. 

2/5 

-1 

4ci|I4| 

'{PT/P^^)maP 

APrlPr^Uin. 

2/5 

+ 1 



( 22 ) 


where Qs = and V^j, = V^/M^ <C 1 is dimensionless. This indicates that if the upcoming 
CMB B-mode polarization experiments would uncover wiggles in the B-mode power spec¬ 
trum, it may offer important information about M^, the string coupling Qs and the shape of 
the dilaton potential. 





FIG. 5: The relation between Prmax/P^^, Mp/{Msy/^) and |F^|. 


In Fig.5, we plot PTmax /with respect to Mp/[Mg^/gl) for different values of |G^|, 
in which ci = 1/8 for the heterotic string. To keep c\ > 0, we require Mp/{Msy/^) < 
-^2/11^01. In addition, the result of Planck [1] showed the scalar perturbation amplitude 
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3 X 10“®, so (5) becomes Mp/{Msy/gl) <C \Mp/H\ ^ 36755 for r = 0.05. In 
performing Fig. 5, all above conditions are satisfied. 

The power spectra plotted in Fig. 3 indicates a'M^ ^ 8 x 10®, 5.4 x 10® and 7.4 x 10® for 
the magenta dotted, green dashed and brown solid curves, respectively, which are consistent 
with Eq. (22) and Fig. 5. Thus Eq. (22) actually works well for the hrst two or three 
oscillations. This provides a direct link between the oscillating amplitude in the primordial 
GWs spectrum and the string parameters, which is depicted in Fig.5. 

In heterotic string theory, the spacetime is 10-dimensional. The leading term of the 
effective action is S'lo = in the string frame, in which 

After compactifying it on a Calabi-Yau manifold with volume Ve, in four-dimension we have 
e.g., [38]. Thus combining it with Fig.5, in which we see Mp /~ 
10"^ for PTmax/Px^ ^ 1-1, we have 


It 


glMl 

M2 


> 1 . 


(23) 


for gs ~ 0.1. This is consistent with the requirement of supergravity approximation, i.e., 
Ve > a'^. 

Here, for simplicity, we required that the background is not affected by the dilaton held. 
However, it is also possible that the dilaton held would ahect the background at certain level. 
In that case, it is required to evaluate the ehect of a' correction on the scalar perturbation 
spectrum. Additionally, the oscillation of primordial GWs power spectrum will also ahect 
the auto-bispectrum of the tensor perturbations and the cross-bispectrum of the primordial 
perturbations, which will be studied in upcoming work. 


IV. CONCLUSION 

The a' correction (1) can be regarded as high-order correction to GR. Due to some 
nonperturbative ehects, the dilaton potential might show itself a local valley. We assume 
that the dilaton potential is hne tuned so that the dilaton held may oscillate rapidly for a 
while around this valley but the inhation background is not ahected. 

We numerically hnd that if the dilaton starts to oscillate at the time of about ~ 60 
e-folds before the end of inhation, the correction (1) may bring unusual oscillations to the 
inhationary GWs spectrum, which might be measurably imprinted in the GMB B-mode 
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polarization. We analytically show that the intensity of oscllating the is determined by the 
string scale Mg and the string coupling gg. 

We show that if r > 0.01, we might observe obvious wiggles around the recombination 
peak in the CMB B-mode power spectrum, which may be well within the detection of the 
current B-mode polarization experiments. This hnding, if supported by observations, would 
yield signihcant insights into the gravity theory beyond GR. Conversely, a null result would 
place tight constraints on corresponding theory. Our work highlights the fact that high- 
precision CMB B-modes polarization experiments might offer us richer information on a 
UV-complete gravity theory than expected. 
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Appendix A: The isocurvature perturbations from cf) 

The curvature perturbation TZ is induced by inflaton, while the perturbation of the spec¬ 
tator held (j) is the isocurvature perturbation. By perturbing 0 with <50, we have 

= j 

- 2ciae-^H^ + 2d^7^(9*50 + TO50) 

+ (Al) 

Thus 50 is instability-free. Actually, we have <C 1, since cia'e~'^H<p 1 and 

0 <C iL, see Eq.(16), otherwise 0 will dominate the background. Thus during inhation the 
curvature perturbation TZ is hardly affected by the spectator scalar held 0. In addition, 
when perturbed, e~'^RQ^ also contributes such as 

~ 7^5050, d^7^(9*7^50, (A2) 

which are higher-order and do not appear in the quadratic action of TZ. This implies that 
as long as inhaton is canonical, TZ is also instability-free. 
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However, after inflation the isocurvature perturbation from 0 possibly converts to the 
curvature perturbation. In this Appendix, we will demonstrate that the isocurvature per¬ 
turbation from 0 is negligible. 

By defining v = a6(j), the perturbation equation of Eq.(12) is given by 


d^v 

dr'^ 



(A3) 


where 2Acia'e~'^H'^{l — e) <C has been neglected and = \ — 

We can set V^f,^ = constant for an analytic estimate [39]. Then the solution can be 
given as 

y (A4) 


When H, the power spectrum of <50 is nearly scale-invariant. However, when ^ H, 

V is imaginary. For convenience, we set v = —iv. Using that hI}\x 1) ~ (a;/2)*^/r(l -|- 
iu) — i(a;/2)*'^cosh(7ri>)r(—zz>)/7r — i(x/2)“*'^r(zz>)/7r, we have \h[^{x l)p for 

D S> 1. Thus the power spectrum of 50 on super horizon scales is given as 


Pfirh — 




I" 


k 


(-), 

\m^J \aH 


(AS) 


which is blue-tilted (~ k^) with a strongly suppressed factor Hence, the isocurvature 
perturbation from 0 is negligible on large scale. 

However, for the quartic potential we used, u varies with time, Eq.(A3) isn’t able to be 
solved analytically. But as long as 10,0 3> H'^, Pst/, is strongly suppressed on the large scale. 
To demonstrate this speculation, we numerically solve Eq.(A3) and plot P 50 in Fig . 6 with the 
same parameters in Fig.l. We can see that the isocurvature perturbation from 0 is negligible 
on large scale. Thus even if after inflation it can be converted to the curvature perturbation, 
its effect to the curvature perturbation from inflaton can be also safely neglected on the 
corresponding scale. 


Appendix B: The 3-|-l decomposition of 

To derive the perturbation action ( 6 ), we should do the 3-1-1 decomposition of the 
Gauss-Bonnet term 

^ This result was first obtained by Kaixi Feng. For more details, see Feng’s upcoming paper. 
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FIG. 6; The power spectrum of Scf) for e = 0.003, = 2.72 x 10 ci = 1/8, a'= 9.6 x 10®, 


Fo/M^ = 6.5 X 10 B = 10^0, 0* = -2.6, ainuiai = -13.4 and 


= -1.5. 


We follow the ADM metric 

NkN’^-N^ Nj 


9^^l' ~ 


, = 


-N' 


-2 


AT* 

W 


Ni 

N‘2 

ij _ 

I 


(Bl) 


where N is the lapse function, W is the shift vector and 7 jj is the spacial metric. It is 
useful to dehne the normal vector of 3-dimensional hypersurface = riQdt/dx^ = (tiq, 0, 0) 
and Using the normalization = —1, one have no = N, which implies 

n^ = {N, 0, 0, 0), n^ = ^), and the 3-dimensional induced metric, which is orthogonal 

to the normal vector, i.e. H^yiny = 0, can be dehned to be H^y = g^y + n^Uy, 

NkN^ N.\ ( 0 0 

Ni 'jij j yO y*-’ 

Thus using the 3-dimensional variables, R^ypa, Rva and R can be written as 


Hpy = 


(B2) 


R 


^Up(7 


= (3) D X K — K K 


p,p-^ ^ ua 


p,CT-^^ Up 


DyK^(j^Tip {^D pKfjp 


D (jKpp)Tly 


-{DpK^p - DyKpp)n^ + {DpKy 

-\-(^DpI\.Q-y D py^Tlp (^Kfy^i^Kp “j- D pQjy H^lKyp “t" Cl i, Cl p^ Tl pTl q- 

~t~ ^p^p. ^n^pp ~t~ ^p^p^^u^a ~l~ ^a(^u ^n^^ua ~t~ ^u^a)^p^p 

Dfjdp HyiKpfj ClpCl(j^TlyTlp^ (®3) 

Ru. = ^^^Rua + KK,, - 2K^,K: - {DJ<: - DJ<)n, - {DpK^^ - D,K)n, 

+ Dpd^ — CjiKpp + ci^cip)nyna- — {D^ydi, — CnKua + Oj^da)-, (B4) 

R = + K^ - - 2{Dpdf^ - H^^C^K^p + d^dp), (B5) 
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where is the covariant derivative with respect to Cn is the Lie derivative in n^- 
direction, = CnH^u/2 is the extrinsic curvature, and Thus after defining 

G ^up DpKjyp DjyKpp, Zpy D(^iy(ip'j py “j- OipQiy, Z Zpy^ WG huve 

jd2 _ 1^2 A D I D 'nLiupa 

-thQB — 

= ^^^Rt,upa - 4 ^^^Rpy 

+K^ + 7{KpyKP''f - 6 ^^^RKpyRP’^ + 2^^^R K'^ + 16 RpyK^'^K!^ 

+4 ^^'’Rpyp^RpPIC^ - UKpyK'^PRp^K^P + WRRpyR’^PR^ - lOR^RpyR^'^ 

-8R ^^^RpyRp^ - 8^^^RpynPRR^pR^^ + SRpyp^RRR^^RP + SG^p^G^ - AGpypG^^P 

+ 16 ^^'^RpyU^GJ^ - 8 Rpyp^GP’^PR + 8 Rpyp^nPnPZ'^^ - 8 ^^^RpyU^RZ + 8 ^^'^RpyZ^’^ 

-4 ^^^RZ + AZRpyRP'^ + SRZpyRp’^ - SZpyR^^R^ - AR^Z (B6) 

Till now, Rqq has been decomposed and presented with only N, Ni and 7 .y-related variables. 

Then, after tedious derivation, the action (6) of the tensor perturbation can be obtained by 
taking ■jij = {t)[e^]ij, where Rj satisfies dihij = 0 and ha = 0 . 

Appendix C: The evolution of </>(«) with different initial value (pinUmi 



FIG. 7: The evolution of (/>(a) with different initial value (piniUah where (pinUiai = —1-3 for red 
dotted curve, (piniUai = —1.4 for green dashed curve, and (pinUiai = —1.5 for brown solid curve, 
respectively. We have set ainuiai = —13.4 for these three cases. 
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